Spectroscopy of ^^F to probe proton-neutron forces close to the drip line 
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A long-lived ,1^ = 4^ isomer, T1/2 ~ 2.2(l)ms, has been discovered at 643.4(1) keV in the 
weakly-bound ^9 F nucleus. It was populated at GANIL in the fragmentation of a '^^S beam. It 
decays by an internal transition to the — if ground state (82(14)%), by /3-decay to ^®Ne, or 
beta-delayed neutron emission to ^^Ne. From the beta-decay studies of the = 1^ and = 4^ 
states, new excited states have been discovered in ^^''^^Ne. Gathering the measured binding energies 
of the J" = 1+ - 4t multiplet in f F, we find that the proton- neutron n0d^^2i'0dg/2 effective force 
used in shell-model calculations should be reduced to properly account for the weak binding of 
^9 F. Microscopic coupled cluster theory calculations using interactions derived from chiral effective 
field theory are in very good agreement with the energy of the low- lying 1^,2^,4^ states in ^®F. 
Including three-body forces and coupling to the continuum effects improve the agreement between 
experiment and theory as compared to the use of two-body forces only. 

PACS numbers: 21.10.Dr, 21.60.Cs, 21.60.De, 23.35. -|-g, 23.20.Lv 



Introduction.- Understanding the boundaries of the 
nuclear landscape and the origin of magic nuclei through- 
out the chart of nuclides are overarching aims and intel- 
lectual challenges in nuclear physics research [ij. These 
are major motivations that drive the developments of 
present and planned rare-isotope facilities. Studying the 
evolution of binding energies for the ground and first few 
excited states in atomic nuclei from the valley of stability 
to the drip line (where the next isotope is unbound with 
respect to the previous one) is essential to achieve these 
endeavours. Understanding these trends and providing 
reliable predictions for nuclei that cannot be accessed 
experimentally require a detailed understanding of the 
properties of the nuclear force [1, ^ . 

In the oxygen isotopes, recent experiments have shown 
that the drip line occurs at the doubly magic ^^Oie [3- 
0], as 25,260 are unbound The role of tensor and 

three-body forces was emphasized in 0, to account 
for the emergence of the = 16 gap at ^^Oie and the 
'early' appearance of the drip line in the O isotopic chain, 
respectively. On the other hand, with the exception of 
^*F [llj l and ^°F which arc unbound, six more neutrons 
can be added in the F isotopic chain before reaching the 
drip line at ^^F22 [l2|. One can therefore speculate that 



the extension of the drip line between the oxygen and 
fluorine, as well as the odd-even binding of the fluorine 
isotopes, arise from a delicate balance between the two- 



body proton-neutron and neutron-neutron interactions, 
the coupling to the continuum [l^l effects and the three 
body forces [13, dl- 

The study of ^^F, which is bound by only 0.80(12) MeV 
[l6| . offers a unique opportunity to investigate several 
aspects of the nuclear force. The ^^F nucleus can be 
modeled using a simplified single-particle (s.p.) descrip- 
tion as a closed ^^O core plus a deeply bound proton in 
the 7r0d5/2 orbital (SV(^^F)~ -15.1(3) MeV All) plus 
an unbound neutron (^^.(^^O)^ 770+io keV 0) in the 
orbital. This simplified picture arises from the 
fact that the first excited state in ^'^O lies at 4.47 MeV 
[3, Q and the nOd^/2 and i'0d;j^2 single particle energies 
are well separated from the other orbitals. The low-lying 
J'' = 15^,2]'',3i',4f states in ^^p thus arise, to a first 
approximation, from the interactions of nuclcons in the 
TTOd^/2 and 1^0^3/2 orbits. 

Present experimental knowledge concerning the mem- 
bers of the J'" = if , 2f , 3f , 4f muhiplet in ^^p is as fol- 
lows. A J'^ = Ij'" assignment has been proposed in [l^l 
for the ground state of ^^F from the observation that its 

2^ states 
The half-life 



beta decay proceeds to the J'^ = Of, J'^ 
and a tentative = Of state in ^^Ne. 
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of F was found to be 10.2±1.4 ms with a P„ value of 
11±4% [Hi. A mass excess AM of 18.680(80) MeV was 
determined for ^^F in [l^ using the timc-of-flight tech- 
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nique. The J'^ = 2^ state was discovered at 657(7) kcV 
jl9| from the fragmentation of ^^'^^Na nuclei. In addi- 
tion a charge-exchange reaction with a ^^Ne beam was 
used in [2^ to study unbound states in ^^F. In this re- 
action, a neutron capture to the 1^(1^/2 orbital and a pro- 
ton removal from the -^(1^/2 (which are both valence or- 
bitals) are likely to occur leading to the J'^ = ~ A'^ 
states. The resonance observed at 271(37) keV above 
the neutron emission threshold [2^ could tentatively be 
attributed to the J'^ = 2>i in ^^F, as it was the only 
state of the J'^ = 1^ — 4]^ which was predicted to be un- 
bound. With the determination of the binding energies 
of the J'^ = 1^ — states, the only missing information 
is the energy of the — 4^ state. In this Letter, we 



demonstrate that the 4^ state is isomeric and decays by 
competing internal transition and /? decay. Its binding 
energy is determined and those of the if — 2^ states are 
re-evaluated. The comparison of the measured binding 
energies of the J'^ = 1]^ — 4f states with two theoretical 
approaches, the nuclear shell model and Coupled Clus- 
ter (CC) theory, provides a stringent test of the nuclear 
forces, where a large proton-to-neutron binding energy 
asymmetry is present. 

Experiment.- The ^^F nuclei were produced through 
the fragmentation of a 77.6 MeV/A 36gi6+ pj-jj^ary beam 
with a mean intensity of 2 jiKe in a 237 mg/cm^ Be tar- 
get. They were selected by the LISE [2l| spectrometer 
at GANIL, in which a wedge-shaped degrader of 1066 
jjjo. was inserted at the intermediate focal plane. The 
produced nuclei were identified from their energy loss in 
a stack of Si detectors and by their time-of-flight with 
respect to the GANIL cyclotron radio frequency. The 
production rate of ^^F was 6 pps with a purity of 22% 
and a momentum acceptance of 2%. Other transmitted 
nuclei, ranked by decreasing order of production, were 
28Ne, 29Na, 27Ne, 249, 22]v^ so^^^ r^j^^^ ^^^^ 

planted in a 1 mm-thick double-sided Si stripped detec- 
tor (DSSSD) composed of 256 pixels (16 strips in the X 
and Y directions) of 3x3 mm^-each located at the final 
focal point of LISE. This detector was used to detect 
the /3-particles in strips j,i±l following the implanta- 
tion of a radioactive nucleus in a given pixel i. With an 
energy threshold of '~80 kcV in the individual strips, a 
/3-efiiciency of 64(2)% was achieved for which was im- 
planted at central depth of the DSSSD. The /3-efficiency 
has been determined from the comparison of the inten- 
sity of a given 7-ray belonging to the decay of ^^F gated 
or not on a /3-ray. Four clover Ge detectors of the EX- 
OGAM array ^ surrounded the DSSSD to detect the 
7-rays, leading to a 7-ray efficiency of 6.5% at 1 MeV. 

The 7-ray spectra obtained up to 2 ms after the im- 
plantation of a radioactive nucleus are shown in Fig. [Ija) . 
In this frame the upper (middle) spectrum is obtained by 
requiring that ^^F (all except ^^F) precedes the detection 
of a 7 ray. A delayed 7-ray transition at 643.4(1) keV is 
clearly observed after the implantation of ^^¥. The bot- 
tom spectrum of Fig. [1] (a) is operated in similar condition 
than the top one, with the additional requirement that 
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FIG. 1: (Color online) (a): 7-ray spectra obtained up to 2 ms 
after the implantation of (upper spectrum), or after the 
implantation of any nucleus except (middle spectrum). 
The bottom spectrum shows the /3-gated 7-rays following the 
implantation of ^''F. (b) Time spectra between implanted ^^F 
and 7-rays, from which half-lives were deduced. The 643.4(1) 
keV and 4+ 2^ (1499.1(4) keV) transitions have the same 
half-life, while the one gated on the 2j 2^ (1672.5(3) keV) 
transition has a larger half- life, (c); /?-gated 7-ray spectrum 
following the implantation of ^^F up to 30 ms. Symbols and 
colors indicate which lines correspond to the /3-decay of the 
1+ (♦, black) and 4+ (■, red) or to the /3 delayed-neutron 
branch (A, blue). The same color codes are used in the decay 
scheme of Fig. [2] Two lines (*, green) could not be placed in 
the decay scheme of ^®F. 



7-rays are detected in coincidence with a j3 transition. 
As the 643.4(1) kcV is not in coincidence with /3 par- 
ticles it must correspond to an internal transition {IT) 
dc-exciting an isomeric state in ^^F, which has a half-life 
of 2.2(1) ms (see Fig. [T{b)). This isomer is likely the 4+ 
state we are searching for. It either decays directly to 
the l"*" ground state, hereby establishing the 4+ state at 
643.4(1) keV. Alternatively, the 643.4(1) keV energy may 
correspond, but with a weak level of confidence, to the 
657(7) kcV state observed in [1^. In this hypothesis, the 
isomerism of the 4+ state would be due to the emission 
of a very low energy 4+ — ;> 2+ transition (up to 10 keV to 
ensure having a long-lived isomer), then followed by the 
2+ 1+ transition. In either case, the excitation energy 
of the 4+ state lies at approximately 650(10) keV. 

The decay of this 4+ state occurs through a competi- 
tion between an internal transition {IT) and /3-dccay to 
two states in ^^Ne. The half-lives corresponding to the 
IT (2.2(1) ms) as weU as to the 1499.1(4) kcV (2.4(2)ms) 
and 1843.4(8) keV (2(l)ms) peaks of Fig. [U^c) are the 
same. These two transitions are seen in mutual coin- 
cidences, as well as with the 2017.6(3) keV 7-ray, pre- 
viously assigned to the 2^ — >■ Of transition in ^^Ne in 
[3. This estabhshes two levels at 3516.7(4) keV and 
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FIG. 2: (Color on line) Decay scheme obtained from the de- 
cays of the 4+ (red) and 1+ states (black) in ^^F to ^^Ne and 
■^^Ne (blue). Shell model predictions obtained with the USDB 
interaction are shown in the right hand side. 



5360.1(9) keV in ^^Ne as shown in Fig. [2 Following 
the Gamow- Teller /3-decay selection rules the 4"*" isomer 
should mainly proceed to the ~ A'^ state in the vibra- 
tor nucleus which we attribute to the 3516.7(4) keV 
state. 

All other observed transitions in Fig. 1(c) from ^^F 
belong to the decay of the 1+ ground state, as their half- 
lives differ significantly from that of the 4+ isomeric state. 
The two 7-ray transitions at 1672.5(3) keV and 1797.1(4) 
keV were found to be in coincidence with the 2017.6(3) 
keV transition, but not in mutual coincidence. This es- 
tablishes two levels at 3690.1(4) keV and 3814.7(5) keV 
which have compatible half-lives of 7.7(2) ms, and 7.8(5) 
ms, respectively. These states presumably belong to the 
two-phonon multiplet of states J'^ — 0^,2^,4^ among 
which the 3516.7(4) keV one was assigned to J'^ = 4^ 
(see above). Using in-beam 7-ra.y spectroscopy from the 
fragmentation of a ^^S beam [2j|, the feeding of the 
3516.7(4) keV level was the largest, that of the 3689.8(4) 
keV state was weaker, while the state at 3814.7(5) keV 
was not fed. As this method mainly produces Yrast 
states, i.e. states having the highest spin value in a 
given excitation energy range, we ascribe J'^ = 2t to 



the state at 3690.1(4) keV, in accordance with [2J|, and 
J-^ = 0+ to the state at 3814.7(5) keV. The fitting of the 
decay half-lives must include the direct 1^ decay of ^^F 
as well as the partial feeding from the 4^ ^ if tran- 
sitions. This leads to a growth at the beginning of the 
time spectrum (Fig. [T] (b) for the 1673 keV 7-ray) which 
depends on the isomeric ratio R and on the internal tran- 
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FIG. 3: (Color online) Calculated and experimental interac- 
tion energies Int(l — 4) in MeV in '^^F. Shell-model calcula- 
tions are shown in the first column using the USDA or USDB 
interactions, while the third column shows results obtained 
with CC calculations. Experimental results are in the center. 
The thickness of the lines corresponds to ilcr error bar. 



sition coefficient IT. These parameters arc furthermore 
constrained by the amount of the 643.4(1) keV 7-rays ob- 
served per implanted ^^F nucleus, leading to R= 42(8)% 
and /r=82(ll)%. 

The /3 feedings derived from the observed 7-ray in- 
tensities are given in Fig. [2j In the /?-delayed neutron 
branch of ^^F to ^'^Ne, some levels observed in [ll[2i,[2l| 
are confirmed, while a new state is proposed at 3114.1(8) 
keV as the 1413.2(7) keV and 1700.9(4) keV 7-rays are 
in coincidence and the summed 7-ray energy is observed 
at 3116(2) keV. A P„ value of 16(4)% (consistent with 
P„=ll(4)% [ll]) is extracted for ^Sp from the observa- 
tion of the 979.7 keV 7-ray in the grand-daughter nu- 
cleus ^^Na whose branching ratio of 18.1(19)% was de- 
termined in pTj . We therefore adopt a mean value of 
P„=13. 50(40)% for ^Sp. The proposed level scheme and 
branching ratios agree relatively well with the shell-model 
calculation shown on the right side of Fig. [2j 

The discovery of this new isomer has an important con- 
sequence on the determination of the atomic mass of the 
26p ground state as well on the interpretation of the one- 
neutron knock-out cross sections from ^^F of Ref. [2^ . It 
is very likely that the measured atomic mass of Ref. [l^ 
corresponds to a mixture of the ground and the isomeric 
states (unknown at that time). As the ^^F nuclei were 
produced in the present work and that of [l6j in similar 
fragmentation reactions involving a large number of re- 
moved nucleons, we can reasonably assume that the ^^F 
isomeric ratio is the same in the two experiments. The 
shift in the ^®F atomic mass as a function of the isomeric 
ratio R amounts to -6.43 keV/%, which for iT!=42(8)% 
yields -270(50) keV. 

Discussion.- The comparison between the experimen- 
tal binding energies of these states can now be made with 
two theoretical approaches, the nuclear shell model and 
CC theory. The experimental (calculated) interactions 
elements arising from the coupling between a pro- 
ton and a neutron, labeled Int(J), are extracted from 



4 



the experimental (calculated) binding energies BE as 

Int(J) =BE(26F)j -BE(26Ffrec). 

In this expression BE(^^F/ree) corresponds to the bind- 
ing energy of the ^^O+lp+ln system, in which the va- 
lence proton in the orbit and the neutron in the 
orbit do not interact. It can be written as 

BE(26Ff,ee) = BE(25F)5/2+ +B^C^0)^,^^ -BIL^'O)^^ . 

Using the relative binding energy of +0.771^q MeV 
between ^^0 and ^^O , the measured atomic masses 
in ^^F and ^^F [l^j and the shift in energy due to the 
isomeric content (see above) it is found that the exper- 
imental value of Int(l) is -1.85(13) MeV. The values of 
Int(2)= -1.19(14) MeV and Int(4)=-1. 21(13) MeV are 
obtained using the J'^ = 2^ and J'^ — 4^ energies of 
657(7)keV and 643.4(1) keV, respectively. A value of 
Int(3)=-0.49(4) MeV is derived from the energy of the 
J'^ = 3^ resonance with respect to the ^^F ground state. 

In the shell-model calculations of Refs. [l^ H^l, the 
two-body matrix elements corresponding to interactions 
in the sd valence space are fitted to reproduce properties 
of known nuclei. Applying these interactions to nuclei not 
included in the global fits (such as bound and unbound 
states in ^^F) implies that shell-model calculations to- 
wards the drip lines can be viewed as predictions. Due 
to the strong coupling to the continuum, and a likely ab- 
sence of many-body correlations not included in the fits, 
these interactions may fail in reproducing properties of 
nuclei like ^^F. Owing to its simple structure, ^^F pro- 
vides a unique possibility to probe the strength of the 
proton-neutron interaction close to the drip line. The 
wave functions of the J'^ = I'l—A'l states are composed 
of mainly (80 — 90%) pure TTQdz/2 ® vQd^/^ component. 
By calculating all states in the J"^ = 1^ — 4^ multiplet, it 
can be seen in Fig. [3]that the J'^ = 1^ state is less bound 
than calculated by about 17% (8%) and that the multi- 
plet of experimental states is compressed by about 25% 
(15%) compared with the USDA (USDB) calculations. 
This points to a weakening of the residual interactions, 
which caused the energy splitting between the members 
of the multiplet. 

We have also performed microscopic CC [111, [32| cal- 
culations for ^^F. This method is particularly suited for 
nuclei with closed (sub-)shells, and their nearest neigh- 
bors. Moreover, CC theory can easily handle nuclei in 
which protons and neutrons have significantly different 
binding energies. To estimate the T:Qd^/2 ^ ^^d^/2 
teraction energy (Int(J)), we use CC theory with singles 
and doubles excitations with perturbative triples correc- 
tions [H,!!^! for the closed-shell nucleus ^'^O, the particle- 
attached CC method for ^^O and ^^ F IssI and the two- 
particle attached formalism for ^^F [36|. We employ in- 
teractions from chiral effective field theory [s^- The ef- 
fects of three-nucleon forces are included as corrections 
to the nucleon-nucleon interaction by integrating one nu- 
cleon in the leading-order chiral three-nucleon force over 



the Fermi sphere with a Fermi momentum fcp in sym- 
metric nuclear matter (ssj . The parameters recently es- 
tablished in the oxygen chain [l^ are adopted in the 
present work. We use a Hartree-Fock basis built from 
A^max = 17 major spherical oscillator shells with the os- 
cillator frequency /lu; = 24 MeV. This is sufficiently large 
to achieve convergence of the calculations for all isotopes 
considered. Using two-body nucleon-nucleon forces we 
get the ground-state energy of ^^F at —173.2 MeV which 
is underbound by ~ 11 MeV compared to experiment. 
However, the relative spectra for the excited states are in 
fair agreement with experiment (see Fig. [3]). In order to 
account for the coupling to the continuum in ^^F, we use a 
real Woods-Saxon basis for the z^lsi/2 and 1^0^3/2 partial 
waves [3^ . The inclusion of continuum effects and three- 
nucleon forces improve the situation, the ground state 
energy is at -177.07 MeV, and the low-lying spectra is 
in very good agreement with experiment. The J'^ = 3"*" 
state in ^^F is a resonance and to compute this state we 
need a Gamow-Hartree-Fock basis [i^ . We are currently 
working on generalizing the two-particle attached CC im- 
plementation to a complex basis. Therefore, the interac- 
tion energy of the J = 3 state is not shown in Fig. O 
Consistently with the shell-model calculations described 
above, a simple picture emerges from the microscopic CC 
calculations: about 85% of the 1+ — 4+ wave hmctions 
are composed of IsOd-shell components, in which config- 
urations consisting of the 7r0c?5/2 and 1^0^3/2 s.p. states 
play a major role. 

Conclusions.- To summarize, a new J'^ = 4]^ isomer 
with a 2.2(1) ms half-life has been discovered at 643.4(1) 
keV. Its isomeric decay to the J'^ = 1^ ground state 
and /3-decay to the J'^ = 4^ state in ^^Ne were ob- 
served. Gathering the /3-decay branches observed from 
the J'^ — 1^ and J'^ = 4^ states, partial level schemes 
of ^^Ne and ^^Ne were obtained. In addition, the ^^F 
nucleus is a benchmark case for studying proton-neutron 
interactions far from stability. The experimental states 
J = 1+ — 4^ arising from the T^d^/2 ® 1^1^3/2 coupling 
in I^F are more compressed than the USDA and USDB 
shell model results. The experimental J'^ — 1^,2^,4^ 
states are less bound as well. These two effects point to 
a dependence of the effective two-body interaction used 
in the shell model as a function of the proton-to-neutron 
binding energy asymmetry. Coupled-cluster calculations 
including three-body forces and coupling to the particle 
continuum are in excellent agreement with experiment 
for the bound low-lying states in ^®F. 
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